PROCEEDINGS 

°F SCIENCE 



O 

(N 

o 
in 



o 

(N 

T— I 

(N 



X: 



BKT phase transitions in strongly coupled 3D Z{N) 
LGT at finite temperature 



O. Borisenko 

BITP, National Academy of Sciences of Ukraine, 03680 Kiev, Ukraine 



E-mail: oleg@bitp.kiev.ua 



V. Chelnokov 

BITP, National Academy of Sciences of Ukraine, 03680 Kiev, Ukraine 



E-mail: vchelnokov@i . ua 



G. Cortese* 

Instituto de Fisica Teorica UAM/CSIC, Cantoblanco, E-28049 Madrid, Spain, Departamento de 
Fisica Teorica, Universidad de Zaragoza, E-50009 Zaragoza, Spain 



E-mail: cortese@unizar.es 



R. Fiore 

Dipartimento di Fisica, Universita della Calabria, and INFN - Gruppo Collegato di Cosenza 
1-87036 Rende, Italy 
E-mail: 



fiore@cs . infn . it 



M. Gravina 

Department of Physics, University of Cyprus, P.O. Box 20357, Nicosia, Cyprus 
E-mail: [gravina @ucy . ac . cy 



A. Papa 

Dipartimento di Fisica, Universita della Calabria, and INFN - Gruppo Collegato di Cosenza 
1-87036 Rende, Italy 



E-mail: papa@cs.infn.it 



I. Surzhikov 

BITP, National Academy of Sciences of Ukraine, 03680 Kiev, Ukraine 



E-mail: i_van_go@inbox . ru 



We investigate, both analytically and numerically, the phase diagram of three-dimensional Z{N) 
lattice gauge theories at finite temperature for N > A. These models, in the strong coupling limit, 
are equivalent to a generalized version of vector Potts models in two dimension, with Polyakov 
loops playing the role of Z{N) spins. It is argued that the effective spin models have two phase 
transitions of infinite order (i.e. BKT). Using a cluster algorithm we confirm this conjecture, 
locate the position of the critical points and extract various critical indices. 

The 30 International Symposium on Lattice Field Theory - Lattice 2012, 
June 24-29, 2012 
Cairns, Australia 



* Speaker. 



© Copyright owned by the atithor(s) under the terms of the Creative Commons Attribution-NonCommercial-ShitreAlike Licence. 



http://pos.sissa.it/ 



BKT phase transitions in strongly coupled 3D Z{N) LGT at finite temperature 



G. Cortese 



1. Introduction and motivations 

After the discovery of the Berezinskii-Kosterhtz-Thouless (BKT) phase transition [[I], ^ |3|], 
almost 40 years ago, this phenomenon still remains an interesting subject. 

It is widely known that this kind of transition occurs in a variety of two-dimensional (2D) 
systems like, for instance, the most elaborated one represented by the 2D XY model. However, 
there are several indications that the BKT phase transitions are also presents in some 3D lattice 
gauge models at finite temperature. Here we study 3D Z{N) lattice gauge theories (LGT) at finite 
temperature in the strong coupling regime. 

While the phase structure of 3D pure Z{N) LGT for N = 2,3 has been the subject of an 
intensive study, much less is known about the finite-temperature deconfinement transition when 
N > 4. On the basis of the Svetitsky-Yaffe conjecture [Q] that connect critical properties of 3D 
Z{N) LGT with the corresponding properties of 2D spin models, we perform this study in order: 

• to clarify the order of the phase transitions that occur (if they are of BKT-type) and then 
to check the prediction for the magnetic critical index rj and the compatibility with the XY 
value for the index v; 

• to confirm the universality with 2D Z{N) vector models and to provide checking-points of 
universality with 3D SU [N] LGT in the strong coupling region. 

1.1 The Model 

We consider a 3D lattice A = xNf with spatial (temporal) extension L (Nt); x = {xq,x\ ,X2) 
where xq G [0,Nt — I] and xi,X2 € [0,L — 1] represent the sites of the lattice and e„ the unit vector 
in the n-th direction. We denote pt (ps) temporal (spatial) plaquettes. It (Is) temporal (spatial) links 
and periodic boundary conditions on gauge fields are imposed in all directions. The conventional 
plaquette angles s{p) is 

s{p) = Sn{x)+S,„{x + en)-Snix + e,„)-S,nix) . (1.1) 

The partition function can be expressed as 

leA s{l)=oJ ih Pt 
where the most general Z(A'^) -invariant Boltzmann weight with — 1 different couplings is 



Q{s) = exp 



2nk 

L /3p(^)cos— s 



(1.3) 



The Wilson action corresponds to the choice j8p(l) = pp, I5p{k) = 0,k = 2,...,N —I. 

To study the 3D Z{N) LGT in the strong couphng limit {p^ = 0) one can map the gauge model 
to a generalized 2D spin Z{N) model with the action 

2 /Ink \ 

^ = I I I cos — {s{x)-s{x + en)) . (1.4) 

X n=\ k=\ V ■'^ / 
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The effective coupling constants are derived from the coupling constant jS, = j3 of the Z{N) 
LGT, using the following equation (the Wilson action is used for the gauge model): 

/3^ = ^lMe.)cos(^) , (1.5) 



where 



p=0 



Bk 



Eexp[,3cos(f^)]c„s(^). (L7) 



In the strongly coupled 3D Z{N) LGT one expects a scenario with thi^ee phases. Therefore 
two phase transitions must separate these thi^ee phases when N > 4 (BKT for N ^ 5): 

• transition from high-temperature to massless phase with T] = 1/4; 

• transition from massless phase to low-temperature phase with the prediction rj =A/N'^. 

1.2 Observables 

Below are listed all the observables used in this work: 

• the absolute value of the complex magnetization: 

Ml = £ exp (^s{x)] = \ML\e'V; (1.8) 

• the real part of the "rotated" magnetization Mr = \Ml\cos{Ny) and the normalized rotated 
magnetization = cos{N\lf) ; 

• the quantity called "population": 

N ( max,-^o.jv_i(n,) 1 \ 
^^ = aF^(, n)' ^^-^^ 

where is number of s{x) equal to /; 

• the related susceptibilities xf^\ Xi^'^\ Xl^ of the real part of the complex magnetization, of 
the rotated magnetization and of the population Si, respectively; 

• the reduced fourth-order Binder cumulant U^'^^ defined as 



j(M) _ , {\Ml\') 

3(|Ml|2)2' 



the cumulant B^^'^^ defined as 



^(M,) _ {\Mr-{Mr)\' 

(|M«-(Mfi)P)2' 



Br' = :::: .::7,;.; ; a-n) 
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• the helicity modulus [Q] 

i:={e)-L^^{s^) , (1.12) 
where e = jj I<,7>^ cos (0,- - dj) and s = I<i7>, sin(0,- - dj) with 0,- = ^s{i). 

We have simulated models with A'^=5,7,9,13 on lattices ranging from L = 128 to L = 2048. 
2. Results 

2.1 Determination of the critical temperatures 

In order to extract the critical indices we need to locate the critical temperatures. Below we 
list the methods used to do that for the first critical coupling jSi'^ 

(a) we locate the positions of the ppl\L) from the peak of the susceptibility xj^^ of on 
various lattice sizes and we find jSi^' by a fit with the following scaling function (with v equal to 
1/2): 

tf^=^iV-^^ ; (2.1) 
(InL + B)- 

(b) we estimate the crossing point of the Binder cumulant u[^^ versus j3 on different lattices 
or, alternatively, we search for the value of jSi'^ which optimizes the overlap of these curves when 
they are plotted against (/3 - j3i^^)(lnL)i/^ (v=l/2); 

(c) we consider the helicity modulus T near the phase transition and define j8pc^(L) as the 
value of P such that T] (j3) = 1/ (2;r/3T) = 1/4 on the various lattices; we then find jSi'^ through the 
function 

valid under the assumption that the phase transition belongs to the XY universality class. 
To determine jSi^^ we use: 

(s) 

(d) the same as the method (a) using instead the susceptibility xl of the population Sl, 

(e) the same as the method (b) using instead simultaneously the Binder cumulant B^^'^^ and 
the order parameter m^,. 

We show in Figs, [l], ^ and |3|the behavior of some of these observables and the method adopted 
to locate the critical couplings. All Tables with the estimations are collected in [^. 

2.2 Determination of the critical indices 

After the determination of the j3f couplings we can extract some critical indices and check 
the hyperscaling relation J = 2j8/v + y/v, where d is the dimension of the system. For the first 
transition, according to the standai^d FSS theory the magnetization \Ml\ at criticality should obey 
the relation \Mi\ ~ L^^l^ for L large enough. Taking into account the possibility of logarithmic 
corrections we use 

\Ml\ = AL-f^/" In' L and x?^ = AL^/^ln'^L , (2.3) 

where the second formula is for the susceptibility Xl^'^' with y/v = 2 — T] (rj magnetic critical 
index). 
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Figure 1: The two plots in the first fine illustrate the method (a), while those on the second line illustrate the 
method (d). 
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Figure 2: The two plots in the first column illustrate the method (b), while the remaining are illustrations of 
the method (e). 
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Figure 3: Behavior of the heUcity modulus T for two models considered, as an illustration of the method 
(c); the solid line represents the curve T = 1 /(InPrj), rj ^ 1/4. 
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Figure 4: Correlation between xi L^^-^ and the Binder cumulant Z?^ in Z(7) with Nr =4- for rj = 0.25 
(up and left) and for 77 = 0.0816 (up and right) for different lattice sizes. MrL^I^ versus m^, in Z(7) with 
Nt = 2 (down and left) and in Z(7) with Nt —4 (down and right) for 77 = 0.0816 for different values of lattice 
size L. 

We apply the same procedure to the second transition with the difference that the fit with the 
scaling laws Eqs. ( ^j| ) is to be applied to data for Mr and respectively. We found, in general, 

a reasonable agreement with the expectations (all the details and tables can be found in [^]). 

The critical exponent rj , for both transitions in these models can be calculated without knowl- 
edge of the critical temperature, building a suitable universal quantity [^]. We show in Fig. ^ some 
results of the behavior of these RG invariant quantities for the Z(7) model with A^, = 2,4. The 
results for the rj index is consistent with the FSS method. 
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3. Summary and Outlook 

We have determined the two critical couphngs ofZ(A'^ = 5,7,9,13) LGT and given estimates 
of the critical indices T] at both transitions. Our findings support for a[\ N ^ 5 the scenario of 
three phases: a disordered phase at high temperatures, a massless or BKT one at intermediate 
temperatures and an ordered phase, occurring at lower and lower temperatures as N increases. This 
matches perfectly with the N limit, i.e. the 3D U (1) LGT (at jS^ = 0), where the ordered phase 
is absent. We have found that the values of the critical index rj at the two transitions ai^e compatible 
with the theoretical expectations. The index v also appears to be compatible with the value 1/2, in 
agreement with RG predictions (see and [||]). 

Considering the determinations of the critical couplings as a function of A'^, we have also 
conjectured the approximate scaling for Pc^''^\n). We found that jSi'^ converges to the XY value 
very fast, like exp(— aA'^^) and jSi^^ diverges like (see [^). 

On the basis of these results we are prompted to conclude that finite-temperature 3D Z{N) 
LGT for A'^ > 4 undergoes two phase transitions of the BKT type and this model also belongs to the 
universality class of the 2D Z{N) spin models, at least in the strong coupling limit jS^ = 0. 

The new numerical techniques tested in these models can be useful to study other models of 
interest like the 3D SU (N) LGT at finite temperature and in the strong coupling region. 

Acknowledgments 

The work of O.B. was supported by the Program of Fundamental Research of the Department 
of Physics and Astronomy of NAS, Ukraine. The work of G.C. and M.G. was supported in part by 
the European Union under ITN STRONGnet (grant PITN-GA- 2009-238353). 

References 

[1] V. L. Berezinsky, Sov. Phys. JETR 32 (1971) 493. 

[2] J. M. Kosterlitz and D. J. Thouless, J. Phys. C. 6 (1973) 1181. 

[3] J. M. Kosterlitz, J. Phys. C. 7 (1974) 1046. 

[4] B. Svetitsky, L. Yaffe, Nucl. Phys. B210 (1982) 423. 

[5] R. Nelson and J.M. Kosterlitz, Phys. Rev. Lett. 39 (1977) 1201. 

[6] O. Borisenko, V. Chelnokov, G. Cortese, R. Fiore, M. Gravina, A. Papa, I. Surzhikov, Phys. Rev. E 86 
(2012) 051131. 

[7] D. Loison, J. Phys.: Condens. Matter 11 (1999) L401. 

[8] O. Borisenko, V. Chelnokov, G. Cortese, R. Fiore, M. Gravina, A. Papa, Phys. Rev. E 85 (2012) 
021114. 



7 



